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Cullis—Radié determinant of a rectangular matrix 
which has a number of identical columns 


ABSTRACT. In this paper we present how identical columns affect the Cullis— 
Radić determinant of an m x n matrix, where m < n. 


1. Introduction. In 1913, Cullis [4] introduced the concept of determinoid 
of a rectangular matrix and it is probably the first published generalization 
of the determinant of a square matrix. Cullis reserved the name determinant 
for square matrices only, but nowadays this term is also used for rectangular 
matrices. 

Since the original Cullis’s definition of the determinant (determinoid) of 
a rectangular matrix is descriptive and requires some nonstandard terminol- 
ogy to be introduced first, we reformulate the idea of Cullis in the following 
way: 


Definition 1.1 (Cullis [4, §3], reformulated). Let A = {a;;} be an m x n 
matrix with m rows and n columns, with elements a;;, where 1 < i < m, 
1< j <n, and let k = min{m,n}. The determinant of A is defined as the 
sum 


det A = 5° (-1)°'P, 
P 
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over all products P = aj,j,...@i,;, of elements of A taken from different 
rows and different columns with 


k 
o(P) =~ ce y 1) $ (i - y i) 
p=1 1<q<p 1<q<p 
tq Stp JaSJp 
where the expression in square brackets is equal to the sum of horizontal and 
vertical steps one must take to pass from the element aij, to the element 
a11 in the matrix obtained from A by removing rows and columns containing 
the elements a;,;, for q < p. 


Independently of Cullis, in 1966, Radić [11] proposed the following defi- 
nition of the determinant, which turns out to be equivalent to the Cullis’s 
definition (the equivalence follows from [4, §30]). 


Definition 1.2 (Radić [11, Definition 1]). Let A = [A1,..., An] bean mxn 


matrix with n columns Aj,...,Ap, and m < n. The determinant of A is 
defined as the sum 
(1.1) det A = So. (Aate det[Aj,,---s Aa 


1<j1<...<jmín 
where r=1+2+...+m. If m >n, then det A = det AT. 


It is worth noting here that there are also other definitions of the deter- 
minant of a rectangular matrix which are not equivalent to the definitions 
of Cullis and Radić, see for example [2, 3, 5, 9, 18, 19, 21]. 

The Cullis-Radić determinant of an m x n matrix, where m < n, and the 
classical determinant of a square matrix have several common properties, 
see [4, §5, §27, §32] and [11, 16], for example: 

(1) The Cullis-Radić determinant of a matrix is a linear function of its 
rows. 

(2) If a matrix A has two identical rows or one of its rows is a linear 
combination of other rows, then the Cullis—Radi¢é determinant of A 
is equal to zero. 

(3) Interchanging any two rows of a matrix changes the sign of its Cullis— 
Radié determinant. 

(4) Adding a linear combination of rows to another row does not change 
the Cullis-Radić determinant. 

(5) The Cullis-Radié¢ determinant can be calculated using the Laplace 
expansion with respect to a row. 

More algebraic properties, which characterize the Cullis-Radié determinant, 
can be found in [1, 4, 6, 8, 12, 13, 14, 16, 20], and some geometric interpre- 
tations are presented in [7, 10, 12, 13, 14, 15, 17, 20]. 

In this paper, we are going to present how identical columns affect the 

Cullis-Radi¢ determinant of an m x n matrix, where m < n. The rest 
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of the paper is organized as follows. In Section 2, we consider the Cullis— 
Radić determinant of matrices with two identical columns. In Section 3 
and Section 4, we present the results for matrices which have an arbitrary 
number of identical adjacent columns and an arbitrary number of identical 
adjacent pairs of columns. Matrices formed by identical adjacent sequences 
of columns are the subject of Section 5. 


2. Two identical columns. First, we introduce some useful notation and 
recall from [6] how interchanging two columns affects the determinant of 
a rectangular matrix. 

Let A = [Aj,...,An] be an m x n matrix, where m < n. We use the 
following notation: 


(i) Aiaj, where i,j € {1,2,...,n}, i Aj, denotes the matrix obtained 
from A by interchanging columns A; and Aj, 

(ii) D;(A) = |A, ie , Åi—1, Aisi, oe , An], where 7 € {1, PAREGA n} and 
n > 1, denotes the matrix obtained from A by removing the col- 
umn Aj, 

(iii) IK (A) = [A1, ...,Aj_-1, K, Ai, Ai+1;, , Anl, where i € {1, PAPE ny}, 
denotes the matrix obtained from A by inserting the column K 
before the i-th column of A. 


Lemma 2.1 ([6, Theorem 2.6]). Let A be an m x (m+1) matriz. Then for 
all i, j € {1,2,..., m +1} such that i # j, we have 


det A + det Ajesj = 0. 


Lemma 2.2 ([6, Theorem 2.7]). Let A =[A1,...,An] be an m x n matriz, 
wherem <n. Then 


det A + det Aij =2 XO (-1)'H1tet-+im det[A 
1Sji<...<JIm<n 
LJF{ji ojm} 


r J (A eam dA Ap Aal 
1<j1 <... <jm <n 
(iEJ, jJ or i¢ J, jEJ) 
J={min{i, j}, max{i j} jo jm} 
card(J)=0 (mod 2) 


ee eee cal 


where r =1+2+...+m and card(X) stands for the cardinality of X. 


Corollary 2.3 ([6, Corollary 2.8]). Let A be an mxn matriz, where m < n. 
Then for every i € {1,2,...,n—1}, we have 


0 ifne{m,m+ 1}, 


ROO SET) Yea ifn > m+2. 
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The following theorem gives a sufficient condition for a matrix to have the 
determinant equal to zero and shows that the determinant of a matrix which 
has two identical columns can be expressed by determinants of matrices 
obtained from the given matrix by removing one or both identical columns. 


Theorem 2.4. Let A = [A),...,A,] be an mxn matriz, where m < n and 
Ay = Ay for some i,j € {1,2,...,n} such thatt < j. Then 
0, ifn €{m,m + 1}, 
det D;(D;(A)), ifn >m+2 andj =i+1, 
det A = jal : 
i 2 E (-1)P-*1 det D:(Dp(A)) 
pait1 
+(—1) i1 det Di(D;(A)), ifn >m+2 andj >i+1. 


Proof. If n =m-+1, then we apply Lemma 2.1 and obtain 
2det A = det A + det Ais; = 0. 
In the case where n > m+ 2 and j = i + 1, by Corollary 2.3, we have 
2det A = det A + det Ain; = 2 det D;(D;(A)). 


Ifn > m+2 and j > i+ 1, then for every k € {0,1,...,j —i— 2} 
Corollary 2.3 yields 


(2.1) det Ii, (D;(A)) = 2det Dj(Di+n41(A)) — det IA%,, , 


From the previous case it follows that 


det I (D,(A)) = det D;(D;(A)). 


(Di(A)). 


Therefore, using (2.1), we obtain 
det A = det IA‘(D;(A)) = 2 det Dj(Di41(A)) — det IA* (Di(A)) 
= 2 det Dj(Dj41(A)) — 2 det Dj(Dj+2(A)) + det 14 (D;(A)) 


j-l 
=2 X. (-1)?-*1 det Di(Dp(A)) + (1) det Di(D;(A)). 
p=i+1 
Remark 2.5. In the proof of Theorem 2.4, instead of (2.1) the following 
relation: 


det I;,(Dj(A)) = 2 det D;(D;j-r-1(4)) — det 1,1 (D;(4)) 


for k € {0,1,...,7 — i — 2} can also be used, yielding 


j—1 
det A=2 $. (-1)-?"! det Dj(Dp(A)) + (—1)7 det D;(D;(A)) 
p=i+1 
forn >m+2andj7>i+1. 
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Example 2.6. Let A = [Aj, A2, A3, A4, A5, Ag] be a 3 x 6 matrix. 
(a) If A2 = As, then 
det A = 2 det D2(D3(A)) — 2 det D2(D4(A)) + det D2(D5(A)) 
= 2det| A1, A4, As, A6] — 2 det[Ai, A3, As, Ae] + det[A1, A3, Aa, A6] 
and 
det A = 2 det Ds5(D4(A)) — 2 det D5(D3(A)) + det Ds(D2(A)) 
= 2det| A1, A2, A3, Ag] — 2 det[ A1, A2, A4, Ag] + det| Aj, A3, A4, A6]. 
(b) If A> = A3, then 
det A = det[Ai, Aa, As, Ae]. 


3. An arbitrary number of identical adjacent columns. In this sec- 
tion, we present two theorems which follow easily from the fact that, while 
calculating the determinant of an m x n matrix, two adjacent identical 
columns can be canceled if n > m + 2 (see Theorem 2.4). 


Theorem 3.1. Let A = [Aj,..., An] be anmxn matrix, wherem < n. Fix 
i € {1,2,...,n} and replace the column A; in the matriz A with k copies of 
A;, where k > 2, obtaining in this way a matrix B of the form 


B= [Anges Ap Ái,- ig Ái, Athan : , An]. 
—— 
k columns 

Then 

det A, if k is odd, 

det B = 4 0, if k is even and m =n, 

det D;(A), if k is even and m <n. 

Theorem 3.2. Let A = [Aj,...,An] be an m x n matriz, where m < n. 


For each i € {1,2,...,n} replace the column A; in the matrix A with ki 
copies of A;, where ki > 1, obtaining in this way a matrix B of the form 


B= A1,- ., A1, A2, ..., A2,- .. 3, An,- -3 Anl: 
— ee — 
kı columns kə columns kn columns 


If p is the number of all odd integers among ky,...,kn and (kj,,...,kj,) is 


the subsequence of all odd numbers taken from the sequence (k1,..., kn), then 
det A, if p=n (all ki’s are odd), 
det B = 4 det[Aj,,...,Aj;,], ifm<p<n, 
0, ifp<m. 


Example 3.3. If A = [Aj, Ag, Ag, A2, A3, A4] is a 2 x 6 matrix and 
P'S [Bi Bo, Bo, BaBa B3, Bi Bs, Bs, Bs, Bs, Be] 
is a 3 x 12 matrix, then 
det A = det[A1, A2, A3, A4] and det B = det| B1, B3, B4, Bel. 
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4. An arbitrary number of identical adjacent pairs of columns. 


Theorem 4.1. Let A = [Aj,...,An] be an m x n matriz, where m < n. 
Fiz i € {1,2,...,m—1} and replace the pair of columns Ai, Aj41 in the 
matrix A with k copies of the pair A;, Ajsi, where k > 2, obtaining in this 
way a matriz By, of the form 


Br = [A1,.--, Ai—1, Ai, Aigi,---, Ai, Aiti, Aita,---, An]. 
k pairs Ay Ai1 
Then 
i Hee ee ifn € {m,m+ 1}, 
k det A — (k — 1) det Dj(Di41(A)), ifn >m+2. 
Proof. First, we prove by induction that for every positive integer k 
(4.2) det By, = (k — 1) det B2 — (k — 2) det A. 
The cases k = 1 and k = 2 are easy. For k > 3, applying Corollary 2.3, we 
have 
det By+ det[Ai,..., Ai—1, Aip Ai, Ai, Ai4i,.--, Ai, Atti, Ai+e,.--; An] 
(4.3) k—1 pairs Aj Aj41 
= 2det By_}. 
Theorem 3.1 yields 
det By + det Bg- = 2 det By_}. 
Assuming that 
det Br_1 = (k — 2) det By — (k — 3) det A, 
det Br_2 = (k — 3) det Bə — (k — 4) det A, 
we obtain 
det By = 2 det By_1 — det Bk—2 = (k — 1) det B2 — (k — 2) det A. 
Considering (4.3) for k = 2, and then applying Theorem 2.4, we have 
2 det A, ifn € {m,m + 1}, 
2det A — det Di(Di+ı(4)), ifn > m+ 2. 
Now, combining (4.2) and (4.4), we obtain (4.1). 
Example 4.2. Let 
A = [Aj, Ag, A3, Aa] and B = [Aj, A2, A1, A2, A1, A2, A3, Aa]. 
(a) If A is a 2 x 4 matrix and B is a 2 x 8 matrix, then 
det B = 3det A — 2 det| A3, A4]. 
(b) If A is a 3 x 4 matrix and B is a 3 x 8 matrix, then 
det B = 3det A. 


(4.4) det Bog = 
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5. An arbitrary number of identical adjacent sequences of columns 
which form the entire matrix. In this section we prove the following 
theorem. 


Theorem 5.1. Let k be a positive integer. If A is anm x n matriz, where 
m <n, then 


k™/? det A, if m is even, 
(m+1)/2 det A. i . : 
det[A,A,..., A] = k det A, he is odd and n is ve 
ss 0, if m,n are odd and k is even, 


k copies of A K(m-1)/2 det A, if m,n, k are odd. 


Proof. The proof follows immediately from Lemma 5.2, Lemma 5.3 and 
Lemma 5.7, which are proved below. 


Before we start with the first lemma, we introduce some useful notation. 
Assume that 


A = Ags sty An] 
is an m x n matrix with n columns Aj,...,An, M < n, and 
B=[|A,A,..., A] = [B1,..., Bk] = [B1,..., Bnk] 
— ES 
k copies of A 

is an m x nk matrix, where B; = A for every i = 1,2,...,k, and By,..., Bnk 
are the columns of B. 

For every sequence a = (aj,...,a%) of nonnegative integers such that 


ae a; = m, define G,(A) to be the set of m x m matrices which have a4 
columns taken from the matrix B,, followed by az columns taken from the 
matrix Bo, etc., followed by az columns taken from the matrix B — in each 
group of a; columns, the columns are arranged in the same order as they 
were arranged in the matrix 6;,i1=1,2,...,n. 

Moreover, for every matrix M = [M),...,Mm] E€ Ga( A) such that M; = 
B;, for each i = 1,2,...,m, define 


cla, M) = N jg 
i=l 


Lemma 5.2. Let A = [Aj,...,An] be an m x n matrix with n columns 
Aj,..., An, where m < n, and let B = [B,,...,B,] be an m x nk matriz, 
where B; = A for each i =1,2,...,k. Then 


det B = 5 `> (—1)7 +M) det M, 


a=(a1,...,€~) MEGa(A) 
ayt+...t+ap=m 
a,€Z, a,>0, t= 1,250.58 


wherer =1+2+4+...+™. 


Proof. The proof follows easily from Definition 1.2. 
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Let n, m, k be positive integers and n = ees € {0,1}. For any 
sequence (a1,..., ap) of nonnegative integers satisfying ze a; = m, define 
Sr an by the following formulas: 


(a1, 


(5.1b) S(m,0) = 1, 
= 0, if mis even 
5.1d S7 =e tl í 
eae Lec!) fi if m is odd, 
= 0 if m is even 
5.1 S? = ý , 
a ELI) en if m is odd, 
v 
T — _1)t(utn) on 
(5.1f) San) 5 ee Seu ay) 
where u > 2, v > 2 are integers and u +v = m, 
(5.1g) Plain sp Ouel 7 Sar, nap)? 
where p < k, a1 £0, ap Æ 0, k > 3, 
(5.1h) ars Ak) a) Sn QR)? 
where p < k, a1 =... = ap = 0, ap+ı #0, k > 3, 
(5.1i) eee = Slana sa nar rto) lasaia? 
where a1 Æ 0, ak £ 0, k > 3. 
Lemma 5.3. Let A = |A1,..., An] be an m x n matrix with n columns 
A1,..., An, where m < n, and let B = |B1,..., Bp] be an m x nk matriz, 
where B; = A for each i = 1,2,...,k. For every sequence of nonnegative 
integers a = (a1,..., ak) such that D ai =m, we have 
(5.2) XO (=1) teM) det M = SŽ det A, 


MEGa(A) 


wherer=1+...+mandn= A 


Proof. The proof is done by induction on k and it is divided into nine steps. 
The first five steps (a)—(e) cover the base cases. The remaining steps (f)—(i) 
complete the induction. 

(a) For k = 1, we have B = A and (5.2) follows immediately from 
Lemma 5.2 and formula (5.1a). 

(b) For k = 2 and a = (m,0), we have Ga(A) = Gimo) (4) = Gim (4) 
and SŽ = Sai = Shin): Therefore, applying case (a), we obtain (5.2). 
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(c) For k = 2 and a = (0,m), for every matrix M = [Mj,...,Mn] € 
G,(A) such that M; = Bj, for each i = 1,2,...,m, we have M € Gim) (A) 


and 
m 


cla, M) = Soi = SG —n)+mn = c((m), M) + mn. 
i=1 


i=1 


Since there is a one-to-one correspondence between matrices in Ga( A) and 
Gim) (A), we have 


5 (pice) det M = air 5 (prt) det M 
MEGa(A) M€Gm)(A) 
= (-1)™" det A. 


Now, applying (5.1c), we obtain (5.2). 
(d) For k = 2 and a = (1,m — 1), every matrix M € G,(A) can be repre- 


sented in the form M = [A,,,...,Ap,,], where (p2,...,Pm) is an increasing 
sequence of integers and p; € {1,2,...,n}, i =1,2,...,m. Moreover, for ev- 
ery sequence (p1, ..., Pm) of integers such that (p2,..., Pm) is increasing and 


pi € {1,2,...,n}, 4 =1,2,...,%, the matrix [Ap,,..., Apm] is an element of 
Ga(A). 

Fix a nondecreasing sequence (p1,...,Pm) of integers from {1,2,...,n} 
and consider the sum 


S(pis---sPm) = (Yt det M 


over all matrices M € G,(A) such that Ap,,...,Ap,, are the columns of M. 
Since any of the columns Ap,,...,Ap,, can be the first column of M, we 
have 


r+pit 5 (pitn) 
S(pi,---;Pm) = (-1) i=2 det|Ay,, Áp2; -< -, Apm] 


m 
a r+pjt > (pi+n) 
iz 


+51) E E EEE ENS. A 
j=2 


+55 Pi . 
= ? (—1) 51 det Agni Apm]; if m is odd and pı <... < Pm, 
0, otherwise. 
Therefore 
A, ifmi 
SHO Gert mT Sono me {OO Em BoM 
MEGa(A) E A 0, if m is even, 


and now (5.2) follows easily from (5.1d). 
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(e) For k = 2 and a = (m—1,1), the proof is similar to the previous case. 
Fix a nondecreasing sequence (p1,...,Pm) of integers from {1,2,...,n} and 
consider the sum 


S(P1,--+5Pm) = XO (1) 1CM) det M 


over all matrices M € Ga(A) such that Ap,,...,Ap,, are the columns of M. 
The calculations 


S(pı, tee , Dm) 
m1 r+ yo pit(pj+n) 
= (-1) i” etl A pian Ai se Maree s Apm: Ap;] 
j=1 m—1 
r+ >) pit(pm+n) 
+(-1) 1 det Ape cece Age Vy Apm] 
rt pin 
—2(-l) = det|Ay,,...,Ap,,], if m is odd and pı <... < pm, 
0, otherwise 
yield 
FE Carte acta = YD Stes. sPm) 
M€Ga(A) pi<..<Pm 
= J(-1)"det A, ifm is odd, 
r= (ieee if m is even. 


Thus, applying (5.le), we obtain (5.2). 
(f) Let k = 2 and a = (u,v), where u > 2, v > 2 are integers and 
u+v =m. In this case we prove (5.2) by induction. Assume that for every 


m x n matrix A and for every & = (u — 1,0), where 0 > 2 is an integer and 
u+ =m <m, we have 


XO (71) M) det M = SË det A, 


MEGz(A) 


where r=1+2+...4+™. 
Notice that 


v+1 
5y Caren) det M = D (E 80--m ; 


MeéG,(A) 1<pi<...<pm<n \j=1 
where 
5j(Piy+++;Pm) = X (=1)7 +M) det M 


is the sum over all matrices M € Ga(A) such that Ap, is the first column of 
M and the columns A,,, where i € {1,2,...,m}\{j}, are the other columns 
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of M. Moreover, for 1 < pı < ... < Pm < n, we have 


reo rit > (ri+n) 
Somia) = do (ED iut det[An,,-.-, Arm] 


1<ry<...<ryu<n 
1<ru4i<...<rm<n 
T1>=Pj 


{72,4 rm }={P1;- Pm {Pi} 


EES De raana bee) 
= (n= EIEN Ata ac tout 


I<ri<...<Tru In 
IKfru+1<...<Tm IN 
T1=Pj 


{ro,..5?m}={p1,--.Pm}\{p;5} 


where ¢([pj,12,---,m], [P1,--+;Pm]) is the number of interchanges of two 
adjacent columns that should be performed to obtain [Ap,,...,Ap,,] from 
[Ap Arz- <, Arm]: Since 


CU pig T2,- --,Tm], [P1; ---, Pml) 
=f pros emh Bri SOT utuk ees yal) 
See AC Sewers Pj, T2,- -3 Tus utj, -3 Tml, [P1 eres 2p 
= (j — Dae ape - , Tu, Tutis- ea heal [Pj+1, - - <, Pml), 


we have 


Zj = = (—1) lr 20 rustut grate) [P51 1--9Pm) 


l<re<...<ruc<n 
sry j<ie<rmen 


{125-0 Tu Tutje Tm} = {P+ Pm} 


The number Z; depends only on the number of elements in {pj+1,.-.,Pm} 


and does not depend on the numbers p1,...,Dm, therefore 
v+1 
So (Hy tetet) det M = det A+ S°(-1) FOV (1) 0-0-0) 7, 
MEGa(A) j=1 


For each j € {1,2,. o define 7; = 1+2 +... + (m— Ey and consider 
a (m — j) x n matrix JO = [AY, S AD with columns Au ) cece AD, 
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We have 


det AW) . (—1)n-0-D) Z, 


ryt S Pi pure 
= Py (=1) mH det [AY 


Pj+1? 7 


AQ re OY) Z, 


1<pj4i<..-<pm<n 


rtd rit >, (ritn) 


= Oe yee ees 


(5.3) I<pj+1<...<Pm íN IKr2<...<Tru In 
I<fru+j <- <rm íN 


{12,. u Tutis m} ={Pj+1 Pm} 
x det[AY),..., AD, AO... , AW) 
2.2 Tu Y Tuj? rA Tr 
LOO 5 (EACAN dot M 
MEG (u—1,v-541) (AM) 
= SCu—10—j41) - det AD), 


Finally, applying (5.1f), we obtain 


v+1 
(-1)' te) det M = det A-S0(-1) ODF sayy 
MEGa(A) j=1 


= det A - SP. 


(g) For an arbitrary integer k > 3 and a = (a1,...,ap,0,...,0), where 
p < k and a; # 0, ap # 0, we have Ga( A) = G(ay,...,.a,)(A) and c(a, M) = 
C((a1,...,@p), M) for every M € G,(A). Therefore, applying (5.1g), we 
obtain 


D yee det M = `> (—1)7tellar-0p), M) det M 
MEGa(A) MEG (ay.,...,ap) (A) 


= gn det A = S$” det A. 
(a1,...,p) a 


(h) Let k > 3 be an integer and a = (aj,...,a%), where p < k, ay =... = 
ap = 0 and ap41 # 0. For every matrix M = [M),..., Mm] E Ga(A) such 
that M; = B;, for each i =1,2,...,m, we have M € Gy (A) and 


Ap415++Ak) 


m m 


cla, M) = ` h = Soi — np) + mnp = c((ap+1, - --, ak), M) + mnp. 


i=1 i=1 
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Since there is a one-to-one correspondence between matrices in Ga(A) 


and G nissan) (4), we have 
`> (—1)7telaM) det M 
MEGa(A) 
= (ow y (—1)7 +a ax) M) det M 
MEG(ap41,-ap) (A) 
SAAT Caen det A. 


Now, applying (5.1h), we get (5.2). 
(i) Let k > 3 be an integer and a = (aj,...,a@%), where a; Æ 0 and ap Æ 0. 
Notice that 


y a aa ea merece 


MEGa(A) 1<pi<...<pm<n 


where 
Spr.: pm) = o(- 1+" det M 


is the sum over all matrices M € Ga(A) such that Ap, Ap,,...,Ap,, are the 
columns of M. Let b = (b1,...,bg—1), where b; = a; for i = 1,2,...,k— 2 
and bk—1 = aķk—-ı + ax, For 1 < py <... < pm < n, we have 


k k 
r+ >> rit). aili—l)n 
S(p- Pm) = X 6) A & det[Ar,,---, Arg] 


Isri<...<ra,<n 
1<1ay41<--<1ay tag <n 
Tete cae a <..<Itm <n 
{P10 ?m}={P1 Pm} 


r+ 5S Pi > b;(i—1)n+azgn 
= (—1) i=l i=2 


Z (—1) (Fort tep-2t1 omh Poi +. +b, att Pm) 


et Agar Aa A phe, 


betbp_9? “*Pby +...4+b, 941? °° 


where the last sum is taken over the indices satisfying the following relations: 
Lr Siete, Sy Le e otas Sy oes 
LS ey hy, Aap <a Sm Sy 


fii era tmy = ADise ses Dias 
1 < poy +... hott dss SD AN, 


{Pb Po ie sees on , Dm} = AP Bice act beesd ty tee Tm}. 


54 A. Makarewicz and P. Pikuta 


Since the sum 
Z= bD (—1) C Fortttp- t1" mh Port.. top atn Pml) 


1< Pp, 4...4b,_9 41 ePm Sn 
{rbi +..+bp 2+1 orm} ={Pb] $...+b, ot] Pm} 
1S1o,4...+bp_9 <L- <T Oy 4. bp 4 <n 
1<1py4...4b, 441 A Sf 


does not depend on the numbers p1, ..., Pm, P1,---;Pm (it depends only on 
the number of elements in {Pbi +...+bp—2+1;-- - ; Pm }), we have 


k n-1 
r+) pit SS bili—l)n 
S(pi, see , Pm) = X (—1) E i= 
x det[A,,,...,A A 


24- Tb +...4bp_ 9? 


x (—1)%"Z, 


Poy t..tbp_gtl?***? Asn] 


where the sum is taken over the indices satisfying the following relations: 
Lory <0), Sn, LS rye Se Toib SM ees 5 
1 < Tbi...+bp 3+1 tan Sy ie £ M, 
{rics Tm} = {propm h 
1< Dort...tbp-2+1 <- < Pm <n, 
{Por +...+br-2+1; <5 Pm} = {ro ah bp oie mI 


Assuming that (5.2) is true for any sequence of k — 1 nonnegative integers 
which sum is equal to m, we get 


` (=F) det M = bD (Apron) det M - (=1)%"Z 
MEGa(A) MEG, (A) 
=det A - S? . (—1)*"Z. 
For every matrix M € G(q,_-1,a,) we have 
c((ak-1,ak), M) = c((ak-1 + ak), M) + apn = c((bk-1), M) + apn. 
Therefore, for every bg—1 X n matrix A, calculations similar to those in (5.3) 
yield 
det A-(—1)%"Z = det A - agian) 
and hence , 
(2 = Sh, 
Finally, applying (5.11), we obtain 
XO (-1)r to) det M = det A - SP - (-1)*"Z 
MEGa(A) 


ak)’ 


= det A- St. S? = det A- S?. 


i (ak-1,a%) 
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The next two lemmas give explicit formulas for Sty is ea 


Lemma 5.4. Letn and m be positive integers, n = ICD? andm = pt+q, 
where p,q are nonnegative integers. 


(a) Ifp#q (mod 2), then 


(b) If both p and q are even, then 


(5.5) Sha = ( p ) =( q ) F (; 3 
2 2 2272 


(c) If both p and q are odd, then 


(5.6) Si = 0: 


(p,q 


Proof. The proof will be done by induction. If p < 2 or q < 2, then (5.4), 
(5.5) and (5.6) follow from (5.1b), (5.1c), (5.1d), (5.1e) and are easy to 
verify. For p > 2 and q > 2, in each of the following cases we apply (5.1f) 
first. 


(aı) If p is odd and q is even, we have 


pt+q-1 q/2 
T 2 2j—1 n)an 2j nyen 
Sho ( p-l pfen TI Ea a TE, 
j=l 


(CECE) 


(a2) If p is even and q is odd, then 


q (q-1)/2 
Ti — ilpHn) an — 2j+1)(p+n) an 
Slp,a) = ED f OETI, = 2. ep nig CEEE 
i=0 j=0 
—1)/2 2 ; = 
E S E (5 3 a _ (a 34 1) (1) 
= co Cone 


= (Fue. 
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(b) If both p and q are even, we obtain 


q/2 q/2 p+q—2 aa 
T = 2j(p+n) an = n 2 
Sy = ea = LY a) 


j=0 j=0 2 


i 
~My ya 


(c) If both p and q are odd, we have 


Ti = 2j(p+n) an 2j7+1)(p+n) on 
Shia De [k-n a ee a E 


Zo Clean \4 cae ^) z 


j=0 


T 

È 
i 
ae 


Remark 5.5. In [4, §33 and §34], Cullis proved the following formula 


k 
X (i—i) 
Qde A= XO (71) detA, Ag], 
1<j1 S... Sjn 
where 1 < m < k < n, A = [A;, .. . , An] is an mxn matrix and the numbers 
rm are defined by 


mM N+ _ p2n+1 _ [P 2n+2 _ 
2m 7 Qom+1 = Gon = (") and Qom+1 = 0, 


where m,n are nonnegative integers and m < n. From Lemma 5.4 it follows 
that 


Qh —m = Seni = as 


Lemma 5.6. Let k, m and n be positive integers, n = -CD and m = 


pA a; be a sum of nonnegative integers. 
(A) If all the numbers ay,...,a,% are even, then 


_ dit. -tak 
(5.7) uis ( 2 al 


ar 
Dores 


(B) If there is exactly one odd number ap among a1,..., ap, then 


= aj+...tap—1 each 
5.8 ote =e a Jne n, 
a PEREN 


(C) If there are at least two odd numbers among a1,...,ap, then 


(5.9) S% a0: 


(41 ,.-,0k 
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Proof. The proof will be done by induction on k. 

For k = 1, we consider only (5.7) and (5.8). They easily follow from 
(5.la). For k = 2, formulas (5.7), (5.8) and (5.9) follow from Lemma 5.4. 
(Note that p in (5.4) and p in (5.8) are generally not equal, but they have 
the same parity if k = 2.) For k > 3, consider the following cases. 

(A) If all the numbers aj,...,a, are even, we apply (5.1i) and obtain 


T T nm 


(aisar) S (ainar wc he) (akid) 


ait. +ak—2tak—1 tak ak—-1 tak ait... +ak 

= 2 2 = ia Y 
~ (a ak—2 Qk—-1t4äk Ak-1 Qk — |a A ie 

eg , 2 2 22 Zee 2 


2.2% 2 


(B) If ap is odd for some p € {1,2,...,k} and a; is even for each i € 
{1,2,...,k}\ {p}, consider two subcases. 
(Bi) If p < k—1, we apply (5.1i) and obtain 


ajt...tazp—1 ak—1t4ak 
E E 2 +1 2 
Sarad a ( a | aae aa? se) 
alae aoa 


Sarai (ae, 
= a a (—1) E K 
(igen) 


(B2) Similarly, if p — k = 0 or p — k = —1, we have 


oe ata tan—t : ak—1+ak—l1 E 
m. = n 2 —k+3)n 
Sfant ~~ rs aa J) (em ag DY 


2 2 2 2 


st tena) Gain 
= Ceri: | 


(C) If there are at least two odd numbers among a1,..., ak, let ap and 
aq, where p < q, p,q E {1,2,...,k}, be the two largest odd numbers in 
{a1,..., ak}. Consider two subcases. 

(C1) If p > 1, we apply (5.1i) k — p times and obtain 


Ti On 
Slan, an) _ E E a S A 


T T 
= Oa de asta A Ce eee Te oa) 


The numbers ap and ap+1 +...+ an are odd, so cae ae = 0 and 
consequently Say wi 0. 


(C2) If p = 1, then S? ka 0 and after applying (5.1i) k — 2 


; (a1,a2+...4 
times, we have 


5 (at joey) = 5 (a1 ,02-+.+04)9 (02,03+..+ax) as Papa ea Cee = 0. 


58 


A. Makarewicz and P. Pikuta 


Lemma 5.7. Let k,m,n be positive integers and n = 


(5.10) 
k2, 
is Alm+))/2. 
TETE = 0, 
Pee y(n) /2, 


aiEZ, a20, i=1,2,...,k 


ae Then 


ifm is even, 

ifm is odd and n is even, 

if m,n are odd and k is even, 
if m,n,k are odd. 


Proof. Denote the left-hand side of (5.10) by L and consider two cases. 
(i) If m is even, we apply Lemma 5.6 (C) and Lemma 5.6 (A), and obtain 


L= 


2 


(a1,..-,%) 
ait+...+ap=m 
ai>0, aj=0 (mod 2), i=1,2,...,k 


Sik cat) 


(a1,-.-,%) 
ait+...+ap=m 
ai>0, aj=0 (mod 2), i=1,2,...,k 


es 


it- Fay mi 
2 ae 
a A =k2, 
Qo De 


(ii) If m is odd, by A we denote the set of all sequences of nonnega- 
tive integers (a1,...,@,) containing exactly one odd number and satisfying 


yy a; =m. By Lemma 5.6 (B), we have 
ait+...taz,—1 
2 


a? ee ae (rears 


p=! (a1,...,a% 
G@p=1 (mod 2) 


Jenene 


Se 


p=l (a1, ap)EA S 27°57? 2 ? 
@p=1 (mod 2) 
k- k-D/2, 


if n is even, 


k 
5 (=1)PH . R-D/2, ifn is odd, 
p=1 


kim+1)/2° if n is even, 
= <0, if n is odd and k is even, 
k®™-1)/2. if n is odd and k is odd. 


Example 5.8. If A = [Aj, Ag, A3] is a 3 x 3 matrix and 


B = [A, A, A] = [A1, A2, A3, A1, A2, A3, A1, A2, A3] 


is a 3 x 9 matrix, then 


det B = 3 det A. 
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